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SUMMARY 


This  report  documents  the  development  of  a  mathematical  model  which  represents 
force  coefficients  acting  on  an  airfoil  in  an  unsteady  environment.  In  addition,  the 
model  is  extended  to  account  for  the  three-dimensional  effects  of  radial  flow  expe¬ 
rienced  by  rotor  blade  sections.  The  methods  developed  in  this  report  are  aimed 
at  obtaining  improved  capability  for  predicting  rotor  blade  section  forco  coeffi¬ 
cients,  particularly  at  or  above  stall  conditions.  The  mathematical  model  has 
been  incorporated  into  the  Government's  Rotorcraft  Flight  Simulation  Program 
C-81  (AGAJ71  version). 

The  mathematical  model  was  correlated  with  two-dimensional  oscillating  airfoil 
test  data,  and  the  Government's  Rotorcraft  Flight  Simulation  Program  C-81 
(AGAJ71  version)  was  correlated  with  full-scale  and  model  rotor  test  data.  The 
results  of  the  correlations  indicate  substantial  improvement  in  rotor  prediction 
capability  at  high  speeds  and  rotor  thrust  coefficients  above  stall. 
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FOREWORD 


This  report  presents  the  results  of  rotorcraft  methods  development  conducted 
under  Contract  DAAJ02-71-C-0045,  Task  1F162204AA4101.  The  work  was  per¬ 
formed  by  The  Vertol  Division  of  The  Boeing  Company  and  was  funded  by  the 
Eustis  Directorate,  U.S.  Army  Air  Mobility  Research  and  Development  Laboratory 
(USAAMRDL).  The  contract  covered  the  period  May  1971  through  July  1972. 

Ronald  E.  Gormont  was  the  Vertol  Project  Engineer;  W.  D.  Vann  and  G.T.  White 
monitored  the  contract  for  the  Eustis  Directorate. 

The  scope  of  the  contract  work  included: 

•  review  of  available  analytical  methods  for  representing  unsteady 
aerodynamics, 

•  selection  of  the  most  appropriate  representation  based  on  correlation 
with  existing  two-dimensional  oscillating  airfoil  test  data  and  require¬ 
ments  of  the  USAAMRDL  C-81  Rotorcraft  Flight  Simulation  digital 
computer  program , 

•  development  of  the  selected  representation  for  incorporation  in  the  C-81 
(AGAJ71  version)  program,  and 

•  development  and  incorporation  of  a  method  for  representing  radial  flow 
effects. 

This  effort  was  conducted  as  part  of  USAAMRDL's  continuing  effort  to  improve 
analytical  prediction  methodologies  for  design  and  evaluation  of  rotorcraft. 

The  author  gratefully  acknowledges  the  assistance  of  F.  D.  Harris,  Manager  of  the 
Boeing  V/STOL  Wind  Tunnel,  and  F.  J.  r’arzanin,  Jr. ,  Senior  Engineer  'n  Vertol 
Dhision's  Rotor  Systems  Group,  for  their  technical  comments,  and  Mr.  T.  Moffa, 
Boeing  Computer  Services,  for  incorporating  the  program  revisions  in  the  C-81 
analysis. 

This  report  is  intended  to  provide  a  theoretical  background  of  the  work  conducted 
under  the  contract,  and  to  sevve  as  a  user's  manual  for  operation  of  the  revised 
Rotorcraft  Flight  Simulation  Program. 
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INTRODUCTION 


The  basis  for  helicopter  rotor  loads  and  performance  analysis  was  formed  in  the 
1920's,  when  Glauert  extended  propeller  strip  theory  to  the  case  of  rotary  wings. 
Although  improvements  were  made  to  Glauert's  analysis,  the  assumption  of  linear 
lift  curve  slopes  and  incompressible  drag  coefficients  remained  until  the  advent  of 
digital  computers.  The  computers  permitted  the  inclusion  of  static  airfoil  force 
coefficients  as  functions  of  angle  of  attack  and  Mach  number. 

Rotor  analysis  methods  became  more  sophisticated  when  nonuniform  downwash  and 
blade  elasticity  were  incorporated.  However,  the  use  of  static  airfoil  force  coeffi¬ 
cients  beyond  stall  angles  of  attack  precluded  proper  analysis  of  the  rotor  when  op¬ 
erating  at  high  lift  levels.  The  cause  was  traced  to  the  omission  of  dynamic  effects 
due  to  rate  of  change  of  angle  of  attack.  This  was  first  analyzed  by  Theodorsen  and 
later  verified  by  tests  of  airfoils  undergoing  unsteady  motions. 

In  1969,  F.  D.  Harris3  demonstrated  significantly  improved  correlation  of  theoretic 
prediction  of  rotor  stall  with  rotor  test  data.  He  used  an  analog  computer  program 
that  incorporated  a  functional  representation  of  oscillating  airfoil  data.  Since  analog 
computers  are  not  well  suited  for  day-to-day  production  use  and  they  cannot  simulate 
the  effects  of  nonuniform  downwash  (more  pronounced  a  over  the  rotor  disc) ,  a  digital 
method  was  developed  to  predict  unsteady  airfoil  coefficients. 

The  Eustis  Directorate,  U.  S.  Army  Air  Mobility  Research  and  Development  Lab¬ 
oratory,  funded  the  development  of  methods  (reported  in  this  document)  to  predict 
rotor  performance  in  the  stall  regimes.  This  work  produced  a  mathematical  repre¬ 
sentation  which  was  incorporated  in  the  Government's  Rotorcraft  Flight  Simulation 
Program  C-81  (AGAJ71  version). 

The  C-81  (AGAJ71  version)  digital  computer  program  is  a  general  flight  simulation 
analysis  for  rotorcraft.  It  may  be  used  to  determine  rotorcraft  trim,  performance, 
stability,  stress  loads,  and  aerodynamic  forces  in  both  steady  flight  and  maneuver 
conditions.  The  program  can  be  used  to  analyze  single-rotor,  tandem-rotor 
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and  side-by-side  rotor  helicopter  configurations  as  well  as  isolated  wind  tunnel 
model  rotors.  The  rotor  blade  simulation  includes  a  fully  coupled  aeroelastic 
modal  analysis.  Prior  to  the  work  described  in  this  report,  the  081  program 
had  a  simplified  unsteady  aerodynamics  representation,  but  it  was  not  sufficiently 
rigorous  to  satisfy  the  requirements  of  general  rotorcraft  configurations.  The 
effort  described  in  this  report  produced  representations  for  unsteady  aerodynamics 
(including  stall  hysteresis)  and  radial  flow  effects. 

The  081  (AGAJ71  version)  digital  computer  program  can  be  used  to  analyze  the 
following  rotorcraft  features: 

•  Configuration  -  Single-rotor,  tandem-rotor,  or  side-by-side  rotor  heli¬ 
copter.  Wind  tunnel  model  rotor  simulation. 

•  Flight  condition  -  Steady  flight  or  time-variant  maneuver. 

•  Phenomena  -  Rotor  and  fuselage  trim  condition,  performance,  stability, 
rotor  airloads,  and  stress  loads. 

The  C-81  (AGAJ71  version)  digital  computer  program  can  be  used  to  analyze  the 
following  rotor  features: 

•  Blade  Elastic  Characteristics  -  Modal  analysis  with  felly  coupled  flap- 
chord -pitch  modes. 

•  Nonuniform  Downwash  Approximation  -  Uniform  downwash  with  empirical 
variation  over  blade  radius  and  azimuth. 

•  Rotor  Blade  Aerodynamics  -  The  program  originally  represented  only 
simplified  unsteady  aerodynamics  and  radial  flow  effects.  As  revised 
herein,  the  program  can  represent  rigorous,  comprehensive,  unsteady 
aerodynamics  and  radial  flow  effects.  The  original  program's  capa¬ 
bility  of  computing  lift,  drag,  and  moment  loads  at  20  radial  stations 
and  at  15-degree  azimuthal  increments  is  unchanged. 

•  Rotor  Configuration  -  Analysis  of  articulated,  teetering,  gimballed  or 
hingeless  rotors. 
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The  discussion  which  follows  presents  a  review  of  available  analytical  methods  for 
representing  unsteady  aerodynamics,  selection  of  the  most  appropriate  represen¬ 
tation,  development  of  the  selected  representation  for  incorporation  in  the  C-81 

(AGAJ71  version)  program,  and  development  of  a  method  for  representing  radial 
flow  effects. 

The  unsteady  aerodynamics  methodology  is  correlated  with  two-dimensional  oscil¬ 
lating  airfoil  data  in  Appendix  I.  Appendix  II  contains  correlations  of  the  revised 

C-81  (AGAJ71)  analysis  with  test  data  for  isolated  rotors  and  complete  rotoreraft 
configurations. 
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SURVEY  OF  UNSTEADY  AERODYNAMICS  METHODS 

The  initial  effort  under  the  contract  consisted  of  a  comprehensive  review  and  eval¬ 
uation  of  available  methods  for  the  prediction  of  aerodynamic  force  and  moment 
coefficients  of  oscillating  airfoils.  The  methods  were  evaluated  according  to  these 
criteria.  First  was  correlation  with  available  test  data  from  airfoils  of  varying 
camber  and  thickness  ratio  {specifically  NACA  0006,  NACA  0012,  Vertol  23010- 
1.  58,  and  Vertol  13006-.7).  Second  was  correlation  with  two-dimensional  oscil¬ 
lating  data  (on  the  airfoils  just  named)  for  the  following  condition^: 

-  Below  static  stall  to  above  dynamic  stall. 

-  With  and  without  yawed  airload  data. 

-  With  and  without  the  spanwise  flow  effects. 

The  third  criterion  was  compatibility  with  the  Government's  Rotor  craft  Flight 
Simulation  Computer  Program.  Three  approaches  offered  potentially  acceptable 
solutions.  They  were  the  methods  of  Ericsson  and' Reding^,  Carta  et  al2,  and 
Harris  et  al3.  The  following  paragraphs  summarize  the  significant  points  of  the 
analyses  and  data  correlations  presented  in  the  three  referenced  methods.  The 
references  do  not  address  each  of  the  above  selection  criteria  in  all  cases;  how¬ 
ever,  based  on  the  available  information,  a  judgement  was  made  concerning  the 
merits  and  deficiencies  of  each  method. 

Ericsson  and  Reding1  describe  a  quasi-steady  theory  in  which  time  history  effects 
are  lumped  into  one  discrete  past-time  event,  and  accelerated  flow  effects  (due  to 
pitch  rate)  are  represented  by  an  equivalent  time  lag.  The  method  yielded  accept¬ 
able  correlation  results.  However,  the  approach  uses  a  semiempirical  graphical 
technique  and  is  not  amenable  to  computer  methods.  In  a  more  recent  report, 
Ericsson  and  Reding4  documented  results  of  a  functional  representation  of  their 
method.  However,  the  functional  approach  is  obtained  at  the  expense  of  accuracy. 

4 

Figure  1  presents  the  correlation  results  from  Ericsson  and  Reding  and  shows 
that  the  graphic  technique  produces  good  correlation  with  the  test  results,  but  that 
the  functional  technique  degrades  accuracy.  The  figure  also  shows  that  the  graphic 
technique  correlates  reasonably  well  with  NACA  0012  pitching  moment  data. 

Carta  et  al2  offered  a  technique  of  generalizing  available  sinusoidal  pitch  oscilla¬ 
tion  data  for  an  NACA  0012  section.  Carta  hypothesized  that  the  airfoil  force 
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FUNCTIONAL  TECHNIQUE  (REFERENCE  4,  FIGURE  40) 
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coefficients  could  be  generalized  as  functions  of  instantaneous  angle  of  attack, 
angular  velocity,  and  angular  acceleration.  In  an  attempt  to  isolate  the  dynamic 
contribution  to  the  force  coefficients,  he  subtracted  the  static  force  coefficients 
from  the  measured  two-dimensional  dynamic  test  data.  The  resulting  incremented 
data,  due  to  dynamic  effects,  were  assembled  in  the  form  of  tabular  cross-plots 
for  use  in  computer  analyses.  Figure  2  depicts  a  comparison  of  Carta's  reconsti¬ 
tuted  dynamic  loops  with  the  measured  test  loops.  Although  this  method  provided 
an  improved  representation  of  blade  element  force  coefficients,  there  are  serious 
drawbacks  which  make  it  unsuitable  for  other  airfoils.  The  table  generated  by 
Carta  is  peculiar  to  the  NACA  0012  airfoil  section  and  also  requires  excessive 
computer  storage  space  which  is  undesirable  in  large  computer  simulations  like 
the  C-81  analysis.  Oscillating  airfoil  tests  for  four  different  sections  (Liiva  and 
Gray6,7)  show  that  the  incremental  force  coefficients  between  static  and  dynamic 
tests  are  dependent  on  the  airfoil  section  tested.  (Figures  11  and  13  of  this  report 
illustrate  the  variation  in  stall  delay  of  four  different  airfoils.)  Thus,  to  utilize 
the  Carta  technique  requires  the  generation  of  a  table  for  the  particular  airfoil  of 
interest  which  is  both  undesirable  and  impractical. 

The  approach  utilized  by  Harris  et  al3,  although  semiempirical,  gives  good  corre¬ 
lation  with  section  tests  data  and  is  amenable  to  computer  methods.  Figure  3 
illustrates  Harris'  correlation  with  the  V23010-1.  58  airfoil,  two-dimensional 
oscillating  data.  The  semiempirical  method  is  an  extension  of  theory  developed 
by  Theodorsen5  and  relies  primarily  on  utilizing  static  2-D  section  force  char¬ 
acteristics  and  a  stall  hysteresis  formulation  based  on  observed  oscillating  airfoil 
test  data.  Harris  developed  functional  equations  for  this  method  and  avoids  the 
computer  storage  requirement  inherent  in  Carta’s  method.  Harris’  method  is  also 
sufficiently  general  to  allow  application  to  general  airfoil  types.  An  additional  fea¬ 
ture  pointed  out  by  Harris3  is  that  utilization  of  the  method  in  an  analog  rotor  anal¬ 
ysis  yielded  greatly  improved  correlation  with  existing  model  rotor  and  full-scale 
rotor  test  data.  Figure  4  shows  the  improved  rotor  correlations  obtained  by  Harris. 
Of  particular  importance  are  the  improved  correlations  at  high  rotor  lift  coefficients. 

Based  on  the  criteria  for  selection,  the  method  of  Harris3  was  chosen  as  offering 
the  best  potential  for  yielding  a  successful  representation  of  unsteady  aerody¬ 
namics  in  the  Government's  Rotorcraft  Flight  Simulation  Computer  Program.  In 
order  to  ensure  that  the  Harris  method  satisfied  all  of  the  selection  criteria,  it 
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Figure  2.  Comparison  of  Oscillating  Airfoil  Data  and 
Carta  Reconstructed  Data. 
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was  correlated  with  two-dimensional  oscillating  data  for  the  four  previously  men¬ 
tioned  airfoils.  Appendix  I  contains  the  results  of  the  extonsive  correlation. 

Since  no  data  is  available  for  yawed  oscillating  section  data,  the  effects  of  radial 
flow  were  assessed  by  correlating  the  C-81  analysis  with  measured  rotor  test 
data.  C-81  was  run  with  and  without  radial  flow  effects  in  order  to  determine  the 
impact  on  correlation.  These  data  are  presented  in  Appendix  II 


DESCRIPTION  AND  DEVELOPMENT  OF  SELECTED 
UNSTEADY  AERODYNAMICS  METHOD 


NCOMPRESSIBLE  POTENTIAL  FLOW  SOLUTION  OF  OSCILLATING  AIRFOILS 

The  flow  environment  encountered  by  a  helicopter  blade  element  changes  rapidly 
as  the  blade  traverses  the  azimuth.  At  azimuthal  positions  (4>)  of  0  and  180  degrees, 
the  velocity  components  are  caused  by  the  angular  rotation  of  the  rotorblade.  At 
^  of  90  and  270  degrees,  the  velocity  is  the  sum  of  (90)  and  the  difference  between 
the  (270)  rotational  and  translational  velocities.  The  variations  around  the  azimuth 
can  result  in  blade  elements  encountering  velocities  near-zero  and  near-sonic  in 
one  rotor  revolution  (which  typically  takes  1/4  second  for  full-scale  rotors).  At 
the  same  time  the  velocity  variation  is  occurring,  the  blade  elements  are  subjected 
to  angle-of-attack  variations  from  zero  or  negative  values  to  values  well  beyond 
the  stall  angle.  Most  classical  rotor  analyses  have  ignored  the  time-variant  aspect 
of  rotorcraft  aerodynamics  and  have  assumed  a  quasi-static  relationship  which  al¬ 
lowed  the  use  of  static  airfoil  force  coefficients  in  evaluating  local  aerodynamic 
forces  acting  on  a  blade  element.  Harris**  and  other  observers  noted  that  the  uti¬ 
lization  of  static  airfoil  characteristics  in  conventional  rotor  analyses  did  not  ac¬ 
count  for  the  demonstrated  ability  of  rotors  to  produce  lift  values  well  beyond  the 
maximum  values  predicted  theoretically.  Some  phenomenon  was  enabling  rotor 
blade  elements  to  generate  more  lift  than  could  be  explained  with  two-dimensional 
static  airfoil  data.  Harris’s  goal  was  to. generate  methodology  which  would 
accurately  predict  rotor  performance  throughout  the  normal  operating  range  and 
well  into  the  rotor  stall  regime.  The  following  development  is  based  on  Harris’s 
work  and  extends  it  to  include  second-order  effects  which  were  ignored  in  the 
original  formulation  due  to  the  limited  capability  of  the  analog  computer  being  used. 

Theodorsen’ s'  development  of  unsteady  aerodynamic  theory  is  the  basis  for  the 
methocNdeveloped  herein.  Theodorsen’s  linear,  incompressible  potential  flow 
theory  yields  the  following  equations  for  a  thin  airfoil  oscillating  in  pitch  and 
heaving  in  translation  about  a  mean  pitch  angle  of  zero  degrees: 
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(1) 


2  ®.  ©.  . 

L  =  pb  (Vir  0  +  ?rh  -  7tba  0  )  +^2  ff p Vb  C  (V0  +  h  +  b  (l/2  -  a)  0  ) 

2  ®  2  ,  ®2  -  ©••  ® 

Mp^  =  -pb  ( n(l/2  -  a)  VbO  +  rrb  (1/8  +  a  )  0  -  ajrbh) 

2  ®  .  .  . 

+  2  pVb  n  (a  +  1/2)  C  (Vo  +  h  +  b  (l/2  -  a)  0  ) 


Figure  5  illustrates  the  terms  and  sign  conventions.  The  following  table  explains 
the  circled  reference  numerals. 

I  REFERENCE  KEY  TO  EQUATIONS  1  AND  2 


Reference 

Numerals 

Term  in 

Lift  Equation 

Moment  Arm  for 
Pitching  Moment 

Source 

(1) 

npb2  VO 

b  (1/2  -  a) 

Noncirculatory,  acting 
at  3/4  chord 

(2) 

jrpb2  h 

ab 

j 

Noncirculatory,  acting 
^  at  l/2  chord 

(3) 

tip  b2  a  0 

ab  J 

2*pVb  C 

(V  0  +h  +  b  (1/2  -  a)0) 


The  pitching  moment  also  includes  a 
t  p  b^  0  /8  term  which  reflects  the 
inertial  effects  for  a  cylinder  of  air 
circumscribing  the  airfoil. 

(a  +  l/2b)  Circulatory,  acting 

at  1/4  chord 


Theodorsen's  equations  were  obtained  by  integrating  pressure  distributions  over 
a  thin  airfoil  (flat  plate)  subjected  to  arbitrary  pitch  and  heaving  motion.  The  pres¬ 
sure  distributions  were  obtained  by  specifying  a  system  of  velocity  potentials 
(sources,  sinks  and  rectilinear  flows)  over  the  airfoil  section.  Equations  (1)  and 
(2)  are  the  resulting  integrated  force  equations  for  section  lift  and  section  pitching 
moment  about  the  quarter  chord.  The  terms  in  these  equations  are  of  two  tyjles: 
circulatory  —  those  which  produce  a  trailed  or  shed  wake  (terms  multiplied  by 
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Theodorsen's  C  function)— and  noncirculatory— those  which  are  due  to  oscillating 
pitch  or  heaving  motion  and  do  not  produce  shed  wake.  By  comparing  the  lift  and 
moment  equations  term  by  term,,  an  analogy  can  be  drawn  for  the  effective  location 
of  lumped  components  of  the  integrated  lift  forces.  Referenced  to  the  circled  num¬ 
erals  in  Equations  (1)  and  (2),  the  source  and  analogous  location  of  the  lumped 
lift  components  are  given  in  the  preceding  table. 


For  further  understanding  of  Equations  (1)  and  (2),  the  reader  is  referred  to 
Bisplinghof  et  al  who  devote  considerable  material  to  this  subject. 

Rearranging  terms  and  introducing  chord  c  and  dimensionless  pitch  axis  location 
PA,  the  following  form  is  obtained: 


L=  (2*)  jV2  -f-  ( y  +  “  c  (PA  -  l/2) 


+  (2*)f  V2  c  C  (0  +-|  +  c  (3/4  -  PA)  j) 


(3) 


M 


PA 


.2,  J_ 

V2 


=  -  (2  7T)  f  v2  ((3/4  -  PA)  —•  +  C2  (1/32  +  (PA  -  l/2)2) 

••  *  , 
(PA  -  1/2) -^r)+  ^TTjjV2  c2  (PA  -  1/4)  c  (o  +~  +  c  (3/4 -PA)^) 


(4) 

In  the  above  equations,  Theodorsen's  term  is  a  complex  quantity  and  a  function 
of  reduced  frequency  k 


C  =  F  +  i  G 


where  k  = 


ccj 

2V 


(5) 

(6) 


and  w  is  the  oscillation  frequency. 

Figure  6  presents  values  for  F  and  G  which  were  determined  by  Theodorsen. 
Physically,  the,  F  term  represents  a  lift  deficiency  term  due  to  oscillation  and 
iG  represents  a  phase  lag  between  the  airfoil  oscillatory  motion  and  the  resulting 
aerodynamic  forces. 
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For  the  trivial  case  of  no  pitch  oscillation  and  a  steady  vertical  translation  velocity 
(corresponding  to  inflow  and  flapping  velocity  for  a  rotor),  F  *=  1.  and  G  *  0  and 
Equations  (3)  and  (4)  reduce  to 


L=(2Jr)|V  c  (0+~> 


(?) 


M 


PA  =  (2*>iv2  °2  (PA -1/4)  (0+±) 


(8) 


The  (  e  +  h/v)  quantity  is  the  quasi-steady  angle  of  attack.  Equations  (7)  and  (8)  are 
the  classical  potential  flow  solution  for  lift  and  moment  on  a  thin  airfoil  section. 

When  the  terms  of  Equations  (3)  and  (4)  are  expanded,  the  C  term  introduces 
imaginary  quantities.  However,  the  following  formulation  is  useful  in  performing 
the  expansion.  Since  Theodorsen's  equations  are  for  deviations  from  a  zero  mean 
angle  of  attack  and  are,  therefore,  vibratory  in  nature. 


let 


0  =  9  sin  ojt 


h  =  h,  sin  cat 

v  CO 
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(9) 


0  „  =  «  0  cos  wt  =  i  ca  0 
v  w  v 

0  =  -w20  sincat  =  -ca20 

V  to}  v 


h  y  (jl»^cos  cot  —  1(0  ll 


•«  2  2 
h  -  -  u  h  sin  cot  =  -co  h 
v  ca  v  s 


(10) 


Substituting 


L  =  (2r)|-V2~(—  +  -  c  (PA  -  l/2)~)+  (2»)-|  V2  c  (F  +  iG) 


+  c  (3/4  -  PA)-^) 


(11) 
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(12) 


+  (2* )  “  V2  c2  (PA  -1/4)  (F  +  i  G)  (  0  +  iffili  +  c  (3/4  -  PA)  — 

*  v  y  V  / 


Expanding  yields 


L  =  <2»)|  -  oCPA-1/2) 


(2rr)-|-  V2c  (f  Oy  +  +  F  c(3/4  -  PA)  f  i  GO  v 


G«o 


r^-Gc  (3/4 -PA)—-) 


(13) 


M. 


PA 


p  2c/  °°v  2  o  « 

=  -<2tt)—  (^(3/4 -PA)  —  +  c  (1/32+  (PA-l/2)2)  -£ 
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ch  \ 
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+  (2tt)  ~  V2  c2  (PA  -  1/4)  (F  0  v  +£!^.ny  +  F  c  (3/4  -  PA)  .l”°.v 
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Referring  to  Equation  (10),  and  noting  that 


Combining  Equations  (15)  with  (13)  and  (14)  yields 
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To  this  point,  the  development  has  considered  oscillation  about  a  mean  angle  of 
zero.  Extending  the  theory  to  include  oscillation  about  £  nonzero  mean,  0  Q,  and 
rearranging  terms  gives 

„  0  2  (o  h  0  \ 

-T  (t  +-?  -MPA-A/2,  -x)  + 


0  /  /  Fh 
,)±v  c(  «o+  (f  -  e  (3/4- PA)  G  f)  V  V 


(2 


(-§•  +  c  (3/4-pA)  4-  ]  °v) 


(18) 


Gh 
_ v 

toV 


mpa  ■  <2’> 


.^4  ( 

2  4  \ 


C*  0  2  °v 

(PA  -3/4)  — -c“  (1/32  +  (PA-1/2)  )  —  + 

V  V 


ch 


(PA -1/2) 


:)+  (2,)  i 


V2c2  (PA -1/4) 


(19) 


(«o+(F-o  (3/4 -PA)  -S^)»v  +  ^  + 


Gh 
_ v 

<oV 


(■g+ 0(3/4 -PA)  i)  iyj 


Theodorsen's  F  and  G  functions  act  only  on  oscillatory  terms  which  produce  a  shed 
wake.  Since  0Q  is  a  steady  term,  let  F  =  1  and  G  =  0  for  the  mean  pitch  term. 


Equations  (18)  anu  (19)  are  the  resultant  formulation  for  the  lift  and  moment  on  an 
airfoil  section  which  is  pitching  and  heaving  about  a  mean  pitch  angle  of  0  0.  The 
equations  as  shown  are  for  incompressible  two-dimensional  flow.  The  next  section 
will  modify  these  expressions  to  the  three-dimensional,  compressible  flow  of  a 
rotary  wing. 
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MODIFICATION  TO  COMPRESSIBLE  FLOW  AND  ROTATING  ENVIRONMENT 


Equations  (18)  and  (19)  discussed  previously  related  section  incompressible  lift 
and  pitching  moment  to  airfoil  pitch  and  heaving  rates  and  accelerations.  In  ad¬ 
dition,  the  equations  contain  an  angular  rate  term  w  which  corresponds  to  the 
oscillation  frequency.  For  the  case  of  a  rotating  wing,  the  dominant  frequency  for 
pitch  and  blade  motion  is  a  l/rev  sinusoidal  variation.  Therefore,  the  angular 
frequency  reference  is  equal  to  the  rotor  angular  speed,  S2  ,  and  Theodorsen's  F 
and  G  functions  will  be  evaluated  for  the  l/rev  frequency.  Equations  (18)  and  (19) 
take  the  form 


0  2  c2  '• 

MyV  - 


•  •  • 

f  0  h 

(_L  +  _L  . 

\  v  v2 


c  (PA -1/2) 


v2, 


(27t)^V2c  ^q  +  (F-c  (3/4 -PA)  ££)  ®y  +  +  (20) 


Fh 


Gh 
_ _ v 

ft  V 


“pA  -  (2.)  {  V2 


+(■%  +  c  (3/4  -  PA)  —J  «' 


4(‘ 


c  0 


(PA  -3/4)  — -p  -c2  (l/32  +  (PA-1/2)2)--  + 


2\  °v 


®  P  _2  2 

(PA-1/2)  —  +  (2*)  y  Vc  (PA-1/4) 
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(21) 


Fh 
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Tv  +(ir  +  c(3/4-pA)-)  °v) 
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where 


F  =  / 

G  =  / 


In  the  development  of  Equations  (1)  and  (2),  Theodorsen  computed  the  potential 

•  •  • 

flow  effects  of  airfoil  velocity  distributions  due  to  pitching  ( 0  mid  C  )  and  heaving 
•  •• 

(h  and  h  ).  For  the  case  of  a  rotating  airfoil,  additional  terms  are  encountered. 

Blade  coning  position  produces  a  term  which  is  derived  in  Figure  7.  The  top  in 
Figure  7  shows  the  velocity  vectors  at  the  leading  and  trailing  edges  of  an  air¬ 
foil  which  are  caused  by  blade  rotation  about  the  shaft.  Utilizing  the  small  angles, 
f-  £  and  «  2*  and  resolving  the  vectors  into  components  which  are  both  perpen¬ 
dicular  to  the  shaft  and  in  the  direction  of  the  span,  we  obtain  the  vectors 
and  <  2*  Now  when  these  vectors  are  transferred  to  the  side  view  in  Figure  7, 
they  resolve  into  components  parallel  to  the  span  and  components  perpendicular  to 
the  blade  element  (Vn  0  and  Vjj  4  2  ^  These  latter  components  give  the 
airfoil  a  sense  of  pitch  motion,  as  shown  in  the  velocity  distribution  diagram  in 
Figure  7 .  The  apparent  pitch  rate  is 

0App  =  ft  sin  0  (22) 


To  illustrate  this  effect,  examine  the  extreme  case  of  a  Hade  coned  to  90° 

(sin  P  =1)  and  having  an  angular  velocity  S2.  The  apparent  pitch  rate  ^App  is 
equal  to  12.  In  the  other  extreme,  the  trivial  case  of  P  =  0  yields  0  App  =  0  since 
the  resultant  velocity  due  to  rotation  is  parallel  to  the  section  chord  line. 


For  small  coning  angles 


-  12  <V  V 


W  =  n  "  =  52  <V 


(23) 


where  0  is  the  mean  coning  angle  and  0  v  is  the  instantaneous  deviation  from 
the  mean. 
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SHAFT 


TRAILING  EDGE 


Figure  7 


Effective  Blade  Element  Pitch  Rate  Due  to 
Blade  Coning. 


Another  term  peculiar  to  the  rotary-wing  problem  is  caused  by  steady  terms  in 
the  normal  velocity  component,  h.  The  velocity  component  normal  to  the  blade 
section  (Figure  8)  is  given  by 

h  =  V,,Tmsin  a  cos  B  -  V__  _  sin  B  cos  a  cos  St'  -  (v  +  Z)  cos  B  (24) 
FLT  s  FLT  s  z 


For  small  flapping  angles 

=  VFLTsin  “s  '  VFLT  ^  C0S  “s  °°s  *  -  \  -  & 


(25) 


For  steady  rectilinear  flight,  the  h  term  is  made  up  of  steady  plus  vibratory  com¬ 
ponents  and  can  be  written 


h  = 


h  +  h 
o  v 


VFLT3in  °S- 


“  VFLT  ^  C0S 


a 


s 


cos  -v 

z 


(V) 


>  (26) 


where  v„  is  the  mean  value  of  induced  velocity  and  v7 ,  is  the  instantaneous 
z(o)  z(v) 

deviation  from  the  mean. 


Again,  since  the  /90  and  hQ  terms  from  Equations  (23)  and  (26)  are  steady  terms, 
Theodorsen’s  F  =  1  and  G  =  0  for  them.  Combining  (23)  and  (26)  with  (20)  ana  (21) 
gives  the  final  form  for  the  rotating  environment. 
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MPA=(27r)  i  V2  ’T  (^(PA-3/4)  (ev+Xl/3)  Y  -(l/32+(PA-l/2)2) 

...  2  ch  \ 
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(c  (3/4-PA)  -f-  +  f)  (0y+Xi^v)j 
Rearranging  the  first  group  of  terms  in  Equation  (28)  for  future  use, 
p  2  c3  /  •  i  \ 

MPA=(27r)  2  V  T  <PA~3/4>  (^V^'V  “(PA-l/2> 


(ev  +  n/sv)  — 

y 
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+  2tt  ~  V*  —  ( - -  -  (PA-3/8)  (fl  +  ftfl)  — ^ 
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4  V  4V2 
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Equations  (27)  and  (29)  apply  to  the  case  of  linear  incompressible  aerodynamics, 
hi  order  to  modify  the  expressions  to  nonlinear  compressible  flow,  the  following 
adjustments  are  required. 


First  an  equivalent  angle  of  attack  may  be  defined  by 


“equ 


0q  +  (F-c  (3/4  -  PA) 


v 


Fh  +  h 
v  o 


•  • 

Gh 

— ~  +  c  (3/4  -  PA)  S2  0  +  (c  (3/4-PA)  F  + 


+  a  0 

V  V 


(30) 


Here  it  is  necessary  to  make  a  distinction  between  the  flow  environment  assumed 
by  Theodorsen  and  the  possible  flow  conditions  experienced  by  a  rotorcraft  blade 
element.  See  Figure  9. 


Theodorsen's  assumption  was  that  the  airfoil  experienced  small  deviations  from 
the  steady-state  flow  condition;  i.e.,  h  <<V.  However,  in  the  possible  rotor  en- 
vironment,  h  may  be  large  compared  to  the  total  resultant  velocity.  This  is  true 
of  the  inboard  blade  segments  and  those  near  the  reverse  flow  area.  In  addition, 
small-angle  assumptions  inherent  in  Theodorsen's  equations  may  be  removed  by 
introducing  the  arc  tangent  function  with  respect  to  the  velocity  component  Uq-, 
resulting  in 

“EQU=  °o+  ff-MSA-PA) 


r  * 

-1  \f . 


Gh 


+  tan 


Fh  +  h  +  — — 
V  V  0  f l 


+  c  (3/4-PA)ft  (3  + 


GV 


(c  (3/4  -  PA)  F  +  )  (0 


/ UT 


lox; 


For  the  linear  case,  2  «r  represents  the  lift  curve  slope  and  the  product  2rr 

is  an  equivalent  Cj.  To  introduce  compressible  nonlinear  effects  merely  requires 
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Figure  9.  Blade  Element  Angle  of  Attack  Assumptions, 


entering  equations  or  table  look-up  data  for  the  airfoil  section  of  interest  and  de¬ 
termining  Cj,  Cd  and  Cm  as  a  function  of  and  Mach  number. 

The  resulting  force  equations  are: 


2  s  Yi 

L=  j  PV2cC1  +  (2n)  —  V2  (0y+n/3)  +  -y  -  (PA-l/2) 
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V  v 


The  first  term  in  the  moment  equation  represents  the  basic  test-measured  pitching 
moment  about  quarter  chord  for  the  static  airfoil  section.  This  component  is  cal¬ 
culated  to  be  zero  when  utilizing  potential  flow  theory,  but  is  finite  for  tests  on 
airfoils.  The  second  term,  calculated  by  potential  theory,  represents  the  con¬ 
tribution  to  pitching  moment  due  only  to  pitch-change-axis  offset  from  the  quarter 
chord  (where  the  circulatory  Cj  acts).  The  third  term  reflects  potential  theory  Cj 
contributions  acting  at  3/4  and  l/2  chord;  the  last  term  reflects  the  inertial 
forces  of  an  air  cylinder  having  a  diameter  equal  to  the  chord. 

Another  limitation  must  be  exercised  in  the  event  of  large  angles  of  attack.  For 
this  case,  only  airload  components  normal  to  the  blade  chord  should  be  added  to 
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the  pitching  moment.  For  large  angles,  a  sin  aBE  component  due  to  drag  is 
introduced,  and  the  Cj  contribution  is  limited  by  cos  agE  resulting  in  the  formu¬ 
lation  of  Equation  (34) . 

Further  simplification  of  these  equations  is  accomplished  by  letting 


C1=  (2*)  T  (6v+nn  + 


(PA-l/2)(«  — 

V  . 


L=  j-  P  V2  c  (0,  +  Oj) 


D  -  T  p  y2  c  °d 


M  =  —  P 
PA  2 


~  P  V2  c2  C  +  ~  P  V2  c2 


Cm+TPV  C  Cd(PA-l/4)sinaBE 


+  Y  P  y2  °2  ^(PA-l/4)  Cj  +  (PA -3/4)  cos  a} 


■n  P  c 


hy-  (PA-3/8)  (  Oy  +  n  0 y)c 


Equation  (34)  represents  the  rotor  environment  formulation  for  two-dimen¬ 
sional  flow.  The  final  requirement  in  the  representation  of  the  rotating  environ¬ 
ment  is  to  account  for  three-dimensional  (radial)  flow.  This  is  described  in  detail 
later  under  RADIAL  FLOW  REPRESENTATION,  However,  suffice  to, say  here 
that  the  yawed  flow  produces  two  effects.  The  first  is  an  increase  in  Cqmax  for 
angles  of  attack  beyond  stall: 


max  a  =  o 
cos  a 
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with  the  limitation  that 


Cl 


max  a 


does  not  exceed  the  linear  lift  curve  slope;  i.  e. , 


C1 

mi  ax  , 
A 


< 


a 

linear 


(36) 


The  second  effect  of  yawed  flow  is 


C,  =  C,  +  C .  (sec  A  -1) 
dA  dA  =  0  dSF 


Incorporating  these  into  Equation  (34)  yields 

12  _  , 
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(38) 
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Equations  (38)  are  the  final  formulations  of  forces  acting  on  a  rotorcraft  blade 
element  when  experiencing  steady,  oscillatory  pitch,  oscillating  heaving,  or  com¬ 
bined  motions.  We  shall  now  describe  the  method  of  modifying  these  expressions 
to  account  for  hysteresis  effects  encountered  when  airfoils  oscillate  into  high 
angles  of  attack  where  stall  occurs. 

STALL  HYSTERESIS  REPRESENTATION 

The  unsteady  aerodynamics  representation  just  discussed  may  be  used  directly  to 
estimate  force  coefficients  for  airfoils  oscillating  in  the  linear  aerodynamics  area 
(i.e.  .below  stall).  For  angles  of  attack  above  stall,  the  theory  is  extended  to 
account  for  hysteresis  effects  which  have  been  observed  in  oscillating  airfoil  tests 
(References  6,  7,  and  11). 
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h  - 


12 

Gross  and  Harris  developed  a  method  of  approximating  the  observed  hysteresis 
loops  by  utilizing  two-dimensional  static  wind  tunnel  data  and  an  empirically 
derived  stall  delay  representation.  The  representation  was  subsequently  used  by 

3 

Harris  et  al  to  obtain  improved  predictions  of  rotor  performance.  The  method 
effectively  modifies  an.  instantaneous  blade  element  angle  of  attack  to  a  reference 
angle  for  use  in  entering  two-dimensional  force  coefficient  data.  The  ref¬ 

erence  angle  of  aUack  is  given  by 


where  the  function  Y  has  been  determined  empirically  from  available  oscillating 
test  data.  Positive  a  reduces  o'  and  delays  stall  and  negative  &  increases 
aREF  anc*  sta^*  The  (sig0  ®)  term  carries  the  sense  of  the  stall  delay.  The  Y 
functions,  which  are  a  function  of  Mach  number,  have  been  determined  for  both 
lift  stall  delay  Y^  and  moment  stall  delay  Y  .  It  is  significant  to  note  that  the 
numerical  values  of  the  two  functions  are  different,  indicating  that  the  influence 
of  a  is  not  the  same  for  both  lift  and  moment  stall  delays.  The  following  develop¬ 
ment  summarizes  the  approach  to  be  used  in  treating  stall  hysteresis  effects. 


Figure  1.0  illustrates  a  typical  trace  for  an  airfoil  in  sinusoidal  pitch  oscillation. 
A;s  noted  ir.  the  figure,  an  airfoil  which  is  set  at  a  mean  pitch  angle  just  below 
static  stall  will  encounter  a  stall  condition  on  the  upstroke  of  its  oscillation.  The 
point  at  which  moment  stall  and  lift  stall  occur  may  be  determined,  and  the  re¬ 
spective  values  of  a  and  a  are  available  from  test  measured  data.  Using  a,  the 
function  ^jc<i/2V|  is  computed  for  the  test  point.  This  dimensionless  parameter 
was  selected  by  Grosc  and  Harris12  because  of  its  similarity  to  the  reduced  fre¬ 
quency  parameter  k.  Test  data  for  a  series  of  points  with  various  mean  angles, 
pitch  frequencies,  and  oscillatory  angles  were  used  to  construct  curves  like  those 
in  Figure  11.  Correlations  of  this  type  confirmed  that  the  selected  parameter  does 
normalize  measured  dynamic  stall  delay  angles  for  a  particular  airfoil  and  Mach 
number.  In  these  figures,  the  values  shown  at  \/\ca/2'V\  =  0  represent  the  stall 

angle  for  the  nonoscillating  or  static  test  conditions.  The  incremental  angle  from 

the  point  to  any  given  test  condition  represents  the  measured  stall  delay  due  to  a 

6  7 

effects.  Data  of  this  type  was  accumulated  by  Liiva  and  Gray  ’  for  four  airfoils: 
the  V2301C-1.  58,  NACA  0012  MOD,  V13006-.7  and  NACA  0006. 
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Figure  10.  Typical  Pitch  Oscillation  Test  Data. 


V23010-1.58  AIRFOIL  V13006-.7  AIRFOIL 

MACH  NO.  =  .4  MACH  NO.  =  .4 


Figure  11.  Typical  Dynamic  Stall 


Analysis  of  the  measured  test  data  revealed  that  some  airfoils  (in  particular  the 
6-percent-thick  sections,  V13006-.7  and  NACA  0006)  could  not  be  adequately 
represented  by  a  single  Y  function.  Consequently,  they  are  represented  by  a  com¬ 
pound  slope  formulation  resulting  in  two  Y  values,  Y^  and  Yg.  For  the  thicker 
airfoils,  a  single  value  of  y  was  determined  to  be  adequate  for  representing  the 
variation  in  stall  delay  with  a.  This  latter  result  is  consistent  with  the  work  of 
Gross  and  Harris  since  they  were  primarily  concerned  with  the  behavior  of  the 
V230110-1. 58  section  in  their  original  work. 

For  the  four  airfoils  studied  under  this  contract,  a  generalized  formulation  of  the 

Y  functions  which  is  independent  of  airfoil  section  was  developed.  The  airfoil 
maximum  thickness  ratio,  t/c  is  the  determining  variable  required  to  compute  Y 
functions  which  agree  with  the  measured  test  values.  Equations  (40)  through  (44) 
summarize  the  stall  hysteresis  formulation  as  a  function  of  Mach  number,  and 
airfoil  maximum  thickness  ratio,  t/c.  Figure  12  is  a  schematic  of  the  formulation; 
linear  interpolation  is  used  between  the  computed  end  points  at  Mj^,  MAX  Y2  and 
at  Mn(Y2  =o),  y2  =  °*  The  value  of  at  which  the  slope  changes  from  Yx  to 

Y  g  is  given  in  Equation  (40),  as  illustrated  in  Figure  11. 


Break 


.  06  +  1.  5  (.  06  -  t/c) 


(40) 
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MACH  NO. 


Figure  12.  Schematic  of  Stall  Hysteresis  Formulation. 
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Figure  13  illustrates  the  correlation  of  the  generalized  stall  delay  formulation  with 
measured  test  data  for  the  airfoils  studied. 


Mow  referring  to  Equation  (39),  can  be  represented  as 


Ref 


(L,  M) 


aBE  “  K1  A 


Dynamic 

Stall..  TtT. 
(L,.  U) 


(45) 


where  AaDynamic  gtali  is  computed  from  Equations  (40)  through  (44)  using  the  incre¬ 
mental  dynamic  stall  representations.  This  formulation  has  the  effect  of  reducing 
angles  of  attack  used  to  enter  data  tables  for  positive  a  ,  and  increasing  angle  of 
attack  for  negative  a .  Hence  the  effect  of  a  (dependent  on  sign)  is  to  adjust  the 
reference  angle  of  attack  either  out  of  stall  (positive  a  )  or  deeper  into  stall 
(negative  a ).  The  dynamic  stall  formulation  is  derived  from  stall  delay  angles 
obtained  during  stall  at  positive  a  (increasing  angle  of  attack).  An  initial  formu¬ 
lation  assumed  a  mirror  image  effect  for  negative  a  \  that  is,  A  a  Stall  Dynamic 
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Figure  13.  Stall  Delay  Functions. 
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had  the  same  sensitivity  for  both  positive  and  negative  a .  Subsequent  correlation 
with  two-dimensional  test  results  indicated  that  this  approach  predicted  premature 
stall  for  negative  a  conditions.  As  a  result,  an  empirical  factor  was  intro¬ 
duced  to  adjust  the  sensitivity  for  a  <  0.  Values  of  K-^  =  1.  for  a  >  0  and  =  .5 
for  d  <  0  give  improved  correlation  with  available  two-dimensional  oscillating  test 
data.  Figure  14  shows  the  effect  on  stall  correlation  due  to  the  application  of  the 
Ki  factor.  =  1.  for  both  positive  and  negative  d  gives  a  premature  stall  pre¬ 


diction.  =  1.  for  a  >  0  and  =  .  5  for  a  <  0  minimizes  the  premature  stall 

prediction.  Airfoil  force  coefficients  obtained  at  aRef  are  identified  as  Cj 
C  ,  and  C .  .  Ref 


a  ne£  for  drag  calculations  is  equal  to  a  pef  used  in  moment  calculations.  Drag 
rise  is  assumed  to  have  the  same  delay  as  moment  stall,  since  moment  stall  in¬ 
dicates  a  center  of  lift  shift  due  to  separated  flow,  which  is  coincident  with  drag 
rise.  This  assumption  is  necessary  since  measured  dynamic  drag  data  are  not 
available. 


The  final  force  coefficients  are  computed  as 


where  the  impact  of  stall  hysteresis  on  lift  adjusts  the  effective  lift  curve  slope 
for  above  s&ll.  Note  that  when  aEE  and  aRe£  are  both  below  stall,  the 
effective  lift  slope  and  are  identical  to  the  linear  valu<  -vhich  would  be  obtained 
without  the  stall  hysteresis  representation.  Only  when  aEE  and/or  aRe{.  arc 
above  static  stall  angles  does  the  representation  affect  values.  Similar  com¬ 
ments  apply  to  the  C  and  C  ,  values  computed. 
r  m  a 
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NORMAL  FORCE  COEFFICIENT  (Cn5 


LEGEND: 


f (Hz)  = 

16.34 

K  = 

.061 

%  - 

.4 

ao  = 

7.33° 

Aa  = 

O 

CO 

CO 

TT 

-  OSCILLATING  DATA,  REF  6 

-  STATIC  DATA,  REF  6 

-  UNSTEADY  AERO  REPRESENTATION 

OF  THIS  REPORT 


K.  =  1.0  FOR  a  >  0 
Kj  =  0.5  FOR  a  <  0 

-  UNSTEADY  AERO  REPRESENTATION 

OF  THIS  REPORT 


ANGLE  OF  ATTACK  (a)  *v  DEGREES 

Figure  14.  Influence, of  Empirical  Factor  on  Correlation  at 
Negative  a  . 
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The  final  formulation  is  derived  by  incorpora  ting  the  s  tall  hysteresis  effects 
(Equation  46)  and  the  radial  flow  effects  (Equations  35,  36,  and  37)  into  the  un¬ 
steady  aerodynamics  representation  (Equation  34).  Prom  Equation  (46),  the 
term  CiRej  j m  0  ^  represents  an  effective,  or  dynamic,  lift  curve 
slope  and  replaces  the  lift  slope/  2n,  derived  from  potential  flow  theory. 


Letting 


aEQU  =  [SGe  Equation  (31)j 


5=  J  [t  -(PA-1/2 )(i'v  +  aev)  ~- 


C .  =  c .  +  C,,  (sec  A  -  1) 

d  dRef  dSF 


the  complete  formulation  is 


L  -  i  f  y2  C  I - - - — - A  [a  „  + 

2  r  1 cos  A  ( -  a0  =0  I  i-EQU  J 

\  L  1 


D  =  ~  p  V2  c  C 


1  2  2 
=  T  p  V  c 


C  +  p  V2  c2  (PA-l/4)  C.  sin  a 


m_  ,  2 

Ref 


d  “  “BE 


.  / 


i„„,  \ 


1  ,,2  2  Ref  \) 

+  2  '  *  <aRefL"  aCj=oy  K° EQU 

+  5  (PA-3/4)  J  cos  a 


hv  -  (PA -3/8)  (  ev  +  a/3v)  c 
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The  unsteady  aerodynamics  methodology  is  correlated  with  two-dimensional  oscil¬ 
lating  airfoil  data  in  Appendix  I.  Appendix  II  contains  correlations  of  the  revised 
C-81  (AGAJ71)  analysis  with  test  data  for  isolated  rotors  and  complete  rotorcraft 
configurations. 
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RADIAL  FLOW  REPRESENTATION 

REVIEW  OF  RADIAL  FLOW  METHODS 

Classical  propeller  and  rotor  analyses  assumed  that  blade  element  airloads  and 
performance  could  be  computed  utilising  the  independence  principle;  i.e.,  aero¬ 
dynamic  forces  on  a  yawed  blade  element  occur  only  in  a  plane  perpendicular  to  the 
span  and>are  proportional  to  the  square  of  the  velocity  component  in  that  plane.  In 
reality,  total  velocity  can  be  substantially  higher  than  the  normal  component,  and 
can  act  at  large  sweep  angles  relative  to  the  span  axis.  Recent  publications  on  this 
subject /indicate  that  accurate  rotor  performance  and  airload  prediction  are  com¬ 
promised  when  the  radial  component  of  velocity  is  ignored. 

There  are  generally  two  approaches  being  pursued  in  an  effort  to  estimate  radial 
flow  .effects  on  rotary-wing  aircraft:  (1)  a  theoretical  approach  to  solving  the  be¬ 
havior  of  airfoil  boundary  layer  in  the  presence  of  radial  flow,  and  (2)  empirical 
interpretations  of  available  test  data  for  yawed  airfoil  sections. 

13 

From  the  theoretical  aspect,  McCroskey  and  Dwyer  have  published  a  compre¬ 
hensive  method  of  analyzing  boundary  layer  growth  over  airfoil  sections  in  a  yawed 
flow  environment.  This  publication  is  the  result  of  a  continuing  study,  starting  with 
basic  theories  and  building  a  solid  theoretical  understanding,  in  an  effort  to  under¬ 
stand  the  impact  of  radial  flow  on  rotor  lift  and  drag.  The  following  is  a  summary 
of  the  significant  features  and  limitations  of  their  analysis: 

•  Features 

Perturbation  Analysis 

Small  crossflow  and  quasi-steady  approximations 
Reduces  independent  variables  from  four  to  two 
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Not  useful  for  separation  characteristics 


Suspect  near  separation 

Gives  relative  significance  of  crossflow  and  unsteady  effects 
Linearized  Analysis 

More  relevant  to  problems  involving  separation 

Gives  qualitative  details  of  flow  field 
Numerical  Analysis 

No  simplifying  approximations 

Treats  flow  field  to  point  of  separation 

Evaluates  significance  of  approximations  of  the  perturbation  and 
linearized  analyses 

Limitations 

Treats  only  laminar  incompressible  flow 
Does  not  treat  oscillating  angle  of  attack 
Neglects  trailing  vorticity 

Treats  separation  results  for  3-D  steady  and  2-D  unsteady  flows,  not  3-D 
unsteady  flows  (real  rotor  environment) 

Complete  potential  flow  (serves  as  boundary  condition  at  edge  of  boundary 
layer)  not  available.  Uses  infinite  blade  with  constant  circulation. 

Applications  to  Rotor  Analyses 

Real  flows  have  transition  to  turbulent  boundary  layer  and  attendant  delay 
in  separation.  Thus,  realistic  representation  of  airfoil  characteristics 
would  require  a  laminar-turbulent  analysis. 


Boundary-layer  analysis  internal  to  production  computer  program  would 
be  restrictive  due  to  time  requirements.  Simplified  representation  (prob¬ 
ably  empirical)  would  be  more  desirable. 

Although  the  work  of  McCroskey  and  Dwyer  represents  a  significant  advancement 
in  understanding  the  mechanism  of  radial  flow  effects  on  rotary-wing  boundary 
layers,  the  present  limitations  of  the  analysis  preclude  use  in  a  rotor  analysis  such 
as  the  Government's  Rotorcraft  Plight  Simulation  Program  C-81  (AGAJ71).  In  par¬ 
ticular,  the  requirements  for  three-dimensional  unsteady  flows  and  a  laminar- 
turbulent  analysis  need  to  be  satisfied  in  order  to  make  such  an  analysis  practical 
and  justifiable.  For  rotor  analysis  applications,  a  more  simplified  or  empirical 
method  is  desirable. 

The  impact  of  empirical  approximations  to  radial  flow  effects  has  been  summarized 
by  Harris^’  Utilizing  simplified  approximations  for  representing  radial  flow 
effects,  Harris  has  been  able  to  demonstrate  marked  improvements  in  rotor  per¬ 
formance  and  airloads  predictions  compared  to  results  predicted  by  quasi-static, 
normal  flow  (independence  principle)  rotor  theories.  Conventional  flow  theories 
assume  that  only  the  velocity  component  normal  to  the  blade  span  axis  need  to  be 
considered  when  computing  airloads.  The  significant  features  of  Harris’  method  are 
that  it  accounts  for  increased  lift  capability,  when  stall  is  delayed  due  to  radial  flow 
on  a  yawed  airfoil,  and  for  increased  skin  friction  drag  (above  values  computed  by 
normal  flow  theories),  which  is  calculated  using  the  resultant  velocity  acting  at  the 
blade  element.  Harris'  method  is  an  empirical  approximation  based  on  observed 
experimental  data  and  is  summarized  in  the  discussion  on  THE  EFFECTS  OF 
RADIAL  FLOW  which  follows.  Due  to  the  simplified  numerical  procedures, 
incorporation  of  an  analysis  of  the  type  developed  by  Harris  can  provide  a  signifi¬ 
cant  improvement  in  predictive  capability  with  little  increase  in  required  computer 
time. 

Harris'  empirical  formulation  for  treating  radial  flow  effects  has  been  selected  as 
the  most  practical  method  available.  The  radial  flow  analysis  which  has  been  in¬ 
cluded  in  the  flight  simulation  program  of  this  contract  has  been  described  earlier 
in  this  report  in  MODIFICATION  TO  COMPRESSIBLE  FLOW  AND  ROTATING 
ENVIRONMENT. 


THE  EFFECTS  OF  RADIAL  FLOW 


The  effect  of  radial  flow  on  blade  element  forces  takes  the  form  of  two  distinct 
phenomena.  Radial  flow  affects  section  drag  level  (in  particular,  skin  friction 
drag)  and  the  orientation  of  the  drag  vector,  and  it  affects  maximum  capability. 
The  following  paragraphs  detail  these  results. 

EFFECT  OF  RADIAL  FLOW  ON  DRAG 

Figure  15  illustrates  the  velocity  diagram  for  a  general  blade  element  of  a  rotor 
blade.  For  simplification,  only  velocity  components  in  the  disc  plane  are  shown. 

Classical  rotor  analyses  assumed  that  the  blade  element  force  vectors  could  be 

assumed  to  act  normal  to  the  plane  span  axis  (i.  e. ,  in  the  direction  of  U^)  and 

2  14  * 

were  proportional  to  Ibp.  However,  Harris  concluded  that  this  approach  was 

optimistic  when  calculating  drag  since  the  airfoil  actually  experienced  a  much 

higher  velocity  acting  at  an  oblique  angle.  He  also  found  that  the  skin  friction  drag 

coefficient  (CdSF)  did  not  decrease  in  the  yawed  plane.  This  is  contrary  to  intuitive 

feel  since  in  the  yawed  plane  the  airfoil  section  appears  to  have  a  lower  thickness 

ratio  t/c,  and  it  has  a  higher  Reynolds  number  due  to  its  increased  chord  and 

velocity.  Both  of  these  phenomena  tend  to  reduce  the  skin  friction  drag  coefficient 

of  conventional  two-dimensional  sections.  In  fact,  some  evidence  indicated  the 

Cdor,  actually  increased  for  yawed  flow.  As  a  first  approximation,  the  conclusion 

was  drawn  that  C^gp  is  independent  of  yaw  angle  and  should  be  evaluated  on  the 

basis  of  chord,  velocity,  t/c,  and  R  normal  to  the  span  axis.  Therefore, 

c 


(53) 


where 


—"ft) 


The  pressure  drag  component  (CdD)  is  unaffected  by  the  assumption  that  it  acts  in 

^  15 

the  blade-normal  plane.  Lippisch  and  Beuschausen  documented  pressure  distri¬ 
butions  for  sections  in  yawed  flows  up  to  A  =  40  degrees  and  concluded  that  the 
chordwise  pressure  distributions  were  essentially  unaltered.  Figure  16  illustrates 
typical  chordwise  pressure  distributions  for  a  9-percent-thick  symmetric  hyper¬ 
bolic  section.  Its  angle-of-attack  range  is  limited,  but  its  Mach-number  range  is 
extensive.  Both  show  good  agreement  when  the  pressure  distributions  for  various 
yaw  angles  are  normalized  by  the  dynamic  pressure  normal  to  the  span.  The  pres¬ 
sure  coefficient  corresponding  to  a  local  Mach  number  of  0  or  1  for  the  given  free 
stream  conditions  is  shown  when  it  falls  within  the  range  of  the  ordinate  scale. 

The  drag  vectors  for  a  blade  element  can  therefore  be  derived  as  follows. 

Skin  Friction  Drag  Components 


AD 


SF 


-JL  C 


c  U 


RESULT 


SF 


TOT 


(54) 


In  the  blade-normal  plane: 


AD, 


SF 


=  JL  C 


NORMAL 


,  c  U  cos  A 
dSF  T0T 


(55) 


but 


UTOT  =  VCosA 


AD  =JLc.  cul  sec  A 
SF  2  dgF  T 


(56) 

(57) 


In  the  blade  radial  direction: 


AD, 


-  JLC ,  cu  sin  A 
SFRADIAL  T  dSF  T0T 


(58) 
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where  sin  A  =  Ur/UtOT 


(59) 


P  C 

ADSFradIAL  =  _2  dSF  CUTOTUR 


Pressure  Drag  Components 


P  -  2 

}  =  -£■  C,  cu 

p  .  2d  T 
NORMAL  p 


RADIAL 


Drag  Summary  Equations 


Combining  the  drag  components  normal  to  the  span  axis: 


“normal  =!c,,2t  °d  ‘eoA+Cd 

\  Oi  p 


4- 0  u2„  (C.  +  C  +  C  (sec  A  -l) 

2  T  \dSF  “p  V 


where 


C-,  +  C.  =  C. 

dSF  dp  d2  -  D 


The  combined  terms  C.  +C,  correspond  to  the  total  C,  value  contained  in  con- 

dSF  dP  d 

ventional  two-dimensional  table  look-up  data,. 

IWmr.  =  JC  u2t  fd  +  (sec  A  -  1)|  (64) 


2-D  SF 


In  the  radial  direction: 


DRADIAL  2  °  UTOT  UR  °d. 


Skin  Friction  Drag  Coefficient 

Since  skin  friction  drag  components  are  normally  not  known  foi  specific  airfoils,  a 
general  formulation  (shown  to  agree  with  a  variety  of  section  data)  is  presented  in 


52 


the  following  paragraph.  This  formulation  may  be  used  on  an  incremental  basis  in 
the  drag  equations  just  presented. 

For  thin  airfoils,  skin  friction  drag  may  be  initially  assumed  as  a  function  of  skin 
friction  on  a  flat  plate : 

C.  =  2  C.  (66) 

dSF  1 

where  Cs  is  the  skin  friction  coefficient  for  one  side  of  a  flat  plate.  Hoerner16  de¬ 
termined  that  the  increase  in  skin  friction  due  to  flow  acceleration  over  an  airfoil 
of  finite  thickness  is 


_ SF  _  -  t 

2C(  c 

The  skin  friction  coefficient  is  obtained  . by  adding  equations  (66)  and  (67) 


CJ  =  2  C  fl  +  2  (68) 

dSF  C/ 

Adding  a  pressure  (or  separation)  drag  component,  CdpR,  equal  to  2  Cf  j 


makes  up  the  total  viscous  profile  drag  ; 


Viscous 

Profile 


=  2c,  (i  +  2i  +  100(i)4) 


The  separation  term  is  not  significant  below  t/cmax  values  of  about  20  percent. 
The  above  equation  is  valid  for  airfoils  with  t/cmax  located  near  30-percent  chord 
position.  For  sections  with  t/ c-max  near  40-  to  DO-percent  chord  (laminar  flow 
airfoils),  Hoerner  obtained  a  similar  expression. 


Results  given  by 


'Viscous 

Profile 


=■•  2  C  [1  +  1.2-  +  60  (— 

/  V  C  \C; 


However,  for  airfoil  sections  used  in  typical  helicopter  applications,  maximum 
thickness  ratio  is  generally  less  than  20  percent,  and  it  is  located  near  30-percent 
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chord.  Therefore,  the  profile  drag  is  dominated  by  the  skin  friction  drag  coeffi¬ 
cient,  and  it  reduces  to 


SF 


(71) 


Hoerner  defines  the  flat  plate  turbulent  skin  friction  as 


C/  = 


.044 


R 


1/6 


for  10  <  R„  <  10 


8 


(72) 


which  gives 


c 


Figure  17  compares  this  expression  (t/c  *  .  1)  with  typical  section  data;  The  range 

g 

of  acceptable  agreement  has  a  lower  limit  of  R  ^  1.7  x  10  ,  where  laminar  flow 

c 

effects  begin  to  dominate. 

EFFECT  OF  RADIAL  FLOW  ON  LIFT 

The  effect  of  radial  flow  on  section  C^  is  manifested  as  an  apparent  increase  in 
Ci  capability  when  referenced  to  velocities  in  the  plane  normal  to  the  span 
axis.  At  angles  below  stall,  there  is  no  observable  effect  on  the  combinations  of 
angle  of  attack  and  C,,;  therefore,  the  slope  of  the  lift  curve  is  unaltered. 

However,  at  high  angles  of  attack,  where  unyawed  static  test  data  would  indicate  a 
stalled  cpndition,  the  yawed  section  exhibits  substantial  increases  in  Cj  capability. 
Figure  18  summarizes  the  findings  of  Carta  et  al2,  Purser  and  Spearman2*,  Critzos 
et  al22,  and  Lizak22  with  respect  to  the  effect  of  yawed  flow  on  section  Cj.  The 
important  effect  is  an  improvement  in  the  Cj^^  attainable  at  high  angles  of  attack, 
where  'normal  two-dimensional  section  data  exhibits  a  stalled  condition.  The  data 
of  Figure  18  is  referenced  to  the  velocity  component  normal  to  the  span  axis  and, 
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Figure  18.  Yawed  Flow  Effects  Delaying  Stall. 
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USER'S  MANUAL 


REVISED  INPUT  REQUIREMENTS 

The  revisions  incorporated  in  the  Government's  Rotorcraft  Flight  Simulation  pro¬ 
gram  as  part  of  this  contract  require  few  changes  in  the  input  requirements  of  the 
program.  (Reference  4  describes  the  basic  input  requirements  of  the  program.) 
Additional  inputs  are  required  only  for  the  unsteady  aerodynamics  and  yawed  flow 
versions  of  the  program. 

The  unsteady  aerodynamics  on/off  switch  has  been  retained  from  the  previous 
version  of  this  program.  It  is  controlled  by  input  location  XMR  (26)  for  the  main/ 
forward/ right  rotor  and  by  location  XTR(26)  for  the  tail/aft/left  rotor.  For  both 
rotors,  a  value  of  location  26  <0.  turns  the  unsteady  aerodynamics  representation 
off,  and'  location  26  >0.  turns  the  representation  on.  Any  group  of  on/off  combi¬ 
nations  may  be  used:  unsteady  aerodynamics  may  be  used  on  no  rotors,  both 
rotors,  forward  rotor  only,  or  aft  rotor  only. 

The  yawed  flow  representation  requires  an  additional  control  input  and  is  designated 
YAWFLO  (N).  N  takes  the  value  l  for  the  main/forward/right  rotor  and  the  value 
2  for  the  tail/aft/left  rotor.  YAWFLO(N)  may  take  four  values  of  significance: 

YAWFLO  (N)  =  0.  no  yaw  flow  effects 

=  1.  yaw  flow  effects  on  skin  friction  drag  only 

=  2.  yaw  flow  effects  on  lift  coefficient  only 

=  3..  yaw  flow  effects  on  lift  and  drag  coefficient 

Again,  any  combination  of  control  inputs  may  be  used  for  the  two  rotors. 

In  the  unsteady  aerodynamics  formulation,  values  are  required  for  angle  of  attack 
at  zero  lift  coefficient.  In  order  to  avoid  time-consuming  table  look-up,  these 
values  have  been  represented  by  a  third-order  curve  fit  as  a  function  of  Mach 
number.  The  representation  is  of  the  form 
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ctc  =Q  -  ALFIN1  (N)  +ALFIN2  (N)  *M^ 

+  ALFIN3  (N)  ALFIN4*M^  (77) 

where  N  =  1  for  main/forward/right  rotor 
N  =  2  for  tail/aft/left  rotor. 

Input  units  are  in  degrees  and  are  converted  internally  to  radians. 

For  skin  friction  drag  coefficients  in  the  yawed  flow  representation  and  stall 
hysteresis  effects  in  the  unsteady  aerodynamics  representation,  airfoil  maximum 
thickness  t/c  ratio  is  required.  This  parameter  is  de  gnated  TDIVC  (N)  and  is 
input  as  a  decimal  fraction  (e.  g. ,  for  a  12-percent-thick  airfoil,  TDIVC  =  0. 12). 

All  these  additional  inputs  are  contained  on  one  input  data  card  per  rotor.  The 
card  is  located  after  the  cards  allocated  to  blade  twist  distribution  and  prior  to 
cards  allocated  for  rotor  elastic  mode  shapes.  For  example,  for  the  main/ 
forward/right  rotor  the  additional  input  follows  XMT  (1-21)  and  precedes  XMRMS 
(1-66),  and  for  the  tail/aft/left  rotor  the  additional  input  follows  XTT  (1-21)  and 
precedes  XTRMS  (1-66).  The  data  card  format  is  consistent  with  other  input  and 
is6F10.0.  The  data  input  order  for  this  card  is:  YAWFL/  ALFIN1  ALFIN2 
ALF1N3  ALFIN4  TDIVC.  Reference  24  gives  the  precise  location  of  these 
data  in  the  input  deck  setup. 

PROGRAM  FLOW  FOR  REVISED  AREAS 

The  sequence  of  calculations  performed  in  the  rotor  analysis  section  has  been 
revised  slightly.  Figure  19  shows  both  the  original  and  revised  flow  paths  of  the 
program.  Minor  revisions  are  required  to  introduce  the  unsteady  aerodynamics 
and  yawed  flow  subroutine  (UNSTED)  into  the  computations  when  these  effects  are 
desired.  For  cases  where  unsteady  zero  and  yawed -flow  are  not  used,  the  flow  is 
unchanged,  except  that  subroutine  CMCALC  has  been  deleted.  Unsteady  aero 
effects  on  pitching  moment  have  been  included  in  the  UNSTED  subroutine  developed 
under  this  contract.  Therefore,  CMCALC  is  no  longer  needed. 
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Figure  19.  Flow  Paths  for  the  C-81  Rotor  Analysis. 


NUMERICAL  STABILITY  AND  EFFECT  OF  UNSTEADY  AERODYNAMICS  ON 
C-81  CONVERGENCE 


Throughout  the  correlation  phase  of  this  contract,  no  numerical  instabilities  were 
encountered.  The  combined  rotor  aeroelastic  solution  and  the  unsteady  aero¬ 
dynamics  formulation  obtained  numerically  stable  solutions  to  the  various  rotor 
flight  conditions  shown  in  Appendix  II  of  this  report.  However,  to  reduce  run  time, 
care  should  be  exercised  in  selecting  initial  estimates  of  flight  constants  used  for 
inputs.  Good  estimates  will  result  in  more  rapid  convergence  of  the  program 
solution  and  a  corresponding  reduction  in  computer  run  time. 

The  use  of  the  unsteady  aerodynamics  option  does  result  in  increased  computer 
running  time  for  some  cases.  This  is  due  to  the  complex  interrelationship  of  the 
nonlinear  unsteady  aerodynamics  forcing  functions  and  rotor  blade  dynamic 
responses.  During  correlation  efforts,  it  was  found  that  the  use  of  the  unsteady 
aerodynamic  option  resulted  in  computer  running  times  which  ranged  from  negligible 

increase  above  the  basic  C-81  running  time  to  as  much  as  a  100%  increase  in  run¬ 
ning  time.  Below  dynamic  stall,  the  use  of  unsteady  aerodynamics  does  not  signif¬ 
icantly  improve  the  prediction  capability,  but  it  does  increase  running  time.  Thus 
it  is  recommended  that  the  unsteady  aerodynamics  option  be  used  only  when 
operating  at  lift  and  speed  conditions  which  produce  rotor  stall. 

Use  of  the  radial  flow  option  results  in  negligible  additional  computer  running  time, 
and  provides  a  more  realistic  representation  of  airfoil  lift  and  skin  friction  drag 

coefficients  than  was  available  in  the  previous  C-81  formulation.  It  is  recommended 
that  the  radial  flow  option  be  used  for  all  computer  cases. 
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APPENDIX  I 

CORRELATION  OF  THEORY  AND  TWO-DIMENSIONAL 

P5CTLLATINS~~AIRP0IL  TE5T"'DATA' - 

Appendix  I  contains  a  detailed  correlation  of  available  oscillating  airfoil  test  data 
and  the  theory \developed  in  this  report.  The  test  data  were  extracted  from 
References  6  and  7  and  include  forced  pitch  oscillation  test  data  for  four  airfoil 
sections:  V2S010-1.  58,  NACA  0012  MOD,  V13006-.7,  and  NACA  0006.  In  addition 
some  limited  data  are  available  for  forced  translational  (heaving)  oscillation  of  the 
V23010-1.  58  and  are  included  in  this  section. 

Airfoils  utilized  in  the  correlation  are  typical  of  sections  currently  used  for  heli¬ 
copter  rotor  applications.  The  V23010-1.  58  and  NACA  0012  type  have  been  used  as 
constant  sections  for  the  entire  span  of  rotor  blades  as  well  as  for  inboard  working 
sections  on  variable  geometry  configurations  (e.g. ,  tapered  chord  rotor  blades). 
The  V13006-.7  and  NACA  0006  have  been  used  as  blade  tip  sections  to  alleviate 
compressibility  effects.  The  tests  of  References  6  and  7  were  conducted  at  Mach 
number,  Reynolds  number,  reduced  frequency  and  mean  angles  of  attack  which 
correspond  to  the  conditions  prevailing  on  the  retreating  blades  of  full-scale  heli¬ 
copters.  Although  the  test  models  had  chords  of  8  and  6.  38  inches,  full-scale 
similarity  was  achieved  by  using  a  variable  density  test  facility.  Testing  was  con¬ 
ducted  in  an  effort  to  compile  a  comprehensive  set  of  data  with  which  theories  for 
dynamic  stall  could  be  developed.  Normal  force  Cn  and  pitching  moment  Cm  co¬ 
efficients  in  two-dimensional  flow  were  determined  by  integrating  measured  differ¬ 
ential  pressures  on  upper  and  lower  surfaces' of  the  sections.  Oscillatory  values  of 
lift  coefficient  Cj  and  drag  coefficient  are  rot  available  for  the  tests. 

The  theory  in  SURVEY  OF  UNSTEADY  AERODYNAMICS  METHODS  (earlier  in 
this  report)  has  been  developed  for  the  more  commor  used  lift,  moment,  and 
drag  coefficients.  In  order  to  compare  theory  with  test  Cn,  the  following  trans¬ 
formation  is  used: 

Cr  =  Ci  cos  a  +  C^  sin  a  (78) 

Computed  and  test  values  of  Cm  compare  directly. 

Figures  20  through  23  compare  pitch  oscillation  test  loops  with  predicted  loops 
(derived  using  the  theory  presented  in  DESCRIPTION  AND  DEVELOPMENT  OF 
SELECTED  UNSTEADY  AERODYNAMICS  METHOD)  for  the  four  airfoils. 
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Figure  20.  Continued  (Sheet  3  of  9) . 
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Figure  21.  Dynamic  Loops  for 
Forced  Pitch  Oscil 
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Figure  21.  Continued  (Sheet  2  of  6) 
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Figure  21.  Continued  (Sheet  5  of  6)  . 
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Figure  22.  Continued  (Sheet  2  of  5) . 
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Figure  22.  Continued  (Sheet  4  of  5) 
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Figure  23.  Dynamic  Iioops  for  the  NACA  0006  Airfoil  in  Forced 
Pitch  Oscillation  (Sheet  1  of  6) . 
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Dynamic  data  are  superimposed  over  the  static  force  coefficients  for  each  partic¬ 
ular  airfoil.  The  identifying  test  conditions  are: 

Mean  Pitch  Angle,  Q!0 

Oscillatory  Pitch  Amplitude,  A  a 

Pitch  Frequency,  f  =  ~  V  (79) 


Reduced  Frequency, 
Mach  Number, 


.  Cjd  Cfff 

2V  =  V 

_  V 

Mn  *  rr 

vSOUND 


J 


Pitch  oscillations  approximated  a  sinusoidal  variation  given  by 


a  =  a0  +  Aa  sin  wt 


(80) 


Figures  20  through  23  span  a  Mach  number  range  of  0. 2  to  0.  6  which  comprises 
the  range >tested.  Mean  angles  of  attack  from  below  static  stall  to  above  dynamic 
stall  are  illustrated  for  each  of  the  four  airfoils  mentioned  previously  at  a  Mach 
number  of  0. 4  and  for  the  V23010-1.  58  section  at  a  Mach  number  of  0.  6,  For  the 
remaining  airfoils  at  Mach  numbers  of  0. 2  and  0.  0,  a  test  point  is  included  for 
forced  pitch  oscillation  about  static  stall  a  which  illustrates  correlation  at  dynamic 
stall.  The  reduced  frequencies  are  predicated  on  a  full-scale  l/rev  variation  in 
pitch  angle  for  a  CH-47C  helicopter  rotor. 

The  overall  correlation  between  test  and  theory  in  Figures  20  to  23  is  good  for 
normal  force  coefficient  Cn  and  acceptable  for  pitching  moment  coefficient  Cm. 

Key  points  of  comparison  are  the  Cn  -  a  combinations  at  which  Cm  is  a  maximum 
with  increasing  a,  and  Cm  ~  a  combinations  at  which  Cm  decreases  rapidly  with 
increasing  a.  Such  conditions  are  indicative  of  dynamic  lift  and  dynamic  moment 
stall  for  the  condition  tested.  In  addition,  the  area  enclosed  within  the  loops  is  a 
measure  of  how  well  the  theoretical  and  test  hysteresis  effects  compare.  For  cer¬ 
tain  cases  the  disparity  in  absolute  levels  for  theory  and  test  loops  is  unexplained. 
The  theory  applies  incremental  adjustments  to  static  data  and  generates  loops 
about  the  static  coefficient  lines.  In  general,  the  dynamic  test  loops  behave  sim¬ 
ilarly;  however,  numerous  test  conditions  have  loops  shifted  away  from  static 
levels.  The  shifts  appear  to  occur  randomly  with  positive  or  negative  displacement 
and  generate  no  perceivable  trends.  Since  there  is  no  theoretical  reasoning  to 
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justify  such  departures,  it  is  believed  that  they  probably  reflect  test  measurement 
inaccuracy  or  zero  shifts  in  test  equipment.  This  discrepancy  has  been  noted  by 
other  investigators.  6>  7 

Figure  24  illustrates  the  effect  of  increasing  the  pitch  oscillation  frequency  beyond 
l/rev  variations.  The  figures  reflect  frequencies  which  correspond  to  2  through 
6/rev  for  a  CH-47B  helicopter  rotor.  Up  to  frequencies  of  4/rev  the  correlation 
is  adequate;  however,  gradual  narrowing  of  the  test  Cn  loop  and  increasingly  large 
negative  Cm  spikes  are  not  predicted  by  theory.  At  the  higher  reduced  frequencies, 
the  correlation  gradually  deteriorates.  The  6/rev  frequency  corresponds  closely 
to  the  first  natural  torsion  frequency  of  typical  rotor  blades,  and  the  above  corre¬ 
lation  indicates  that  when  a  blade  is  fluttering  at  its  natural  frequency  with  large 
amplitude  the  airloads  prediction  using  the  developed  theory  would  be  less  than 
ideal.  However,  of  greater  significance  are  the  l/rev  variations  which  dominate 
the  rotor  during  normal  operation  and  which  initially  excite  a  stall  flutter  condi¬ 
tion.  It  is  this  initial  excitation  which  is  of  interest  since  it  is  the  rotor  condition 
to  be  avoided  in  design  and  analysis.  Figures  20  through  23  previously  demon¬ 
strated  that  the  developed  theory  gives  good  correlation  for  the  dominant  l/rev 
pitch  oscillations.  This  aspect  is  explored  further  in  Appendix  n  where  complete 
rotor  configurations  are  correlated  with  the  C-81  theory. 

Blade  spanwise  bending  and  flapping  motion  produce  oscillating  heaving  velocities 
which  act  on  rotor  blade  elements.  Figure  25  is  representative  test  data  for  trans¬ 
lational  oscillation.  For  these  figures  the  identifying  test  condtion,  Ah,  is  the 
oscillating  translation  position  and  is  measured  in  semichords.  The  translation 
oscillations  approximate  sinusoidal  valuations  of  the  type 

h  =  Ah  sin  wt  (81) 

The  theory  and  test  data  compare  favorably,  indicating  that  the  theory  adequately 
treats  velocity  components  due  to  translational  motion. 
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Figure  24.  Continued  (Sheet  2  of  5) 
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Figure  24.  Continued  (Sheet  3  of  5) 
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Figure  25.  Continued  (Sheet  3  of  4). 
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APPENDIX  n 

CORRELATION  OF  THEORY  WITH  ISOLATED  ROTOR 
AND  HELICOPTER  TEST  DATA 


OBJECTIVE 

The  methods  developed  in  this  report  are  aimed  at  obtaining  improved  capability 
for  predicting  airfoil  force  coefficients,  particularly  at  or  above  stall  conditions. 
The  unsteady  aerodynamics  and  radial  flow  methodologies,  presented  herein, 
serve  to  improve  the  prediction  of  rotor  force  coefficients  at  high  rotor  lift  levels 
and  high  forward  speeds. 

Rotor  operation  at  high  advance  ratios  produces  considerable  radial  flow  along  the 
rotor  blade  span.  The  three-dimensional  effects  due  to  radial  flow  are  important 
in  relieving  lift  stall  and  in  increasing  skin  friction  dragiover  the  airfoil  sections. 
In  addition,  high  rotor  lift  and  high  advance  ratios  necessitate  high  lift  coefficients 
on  the  low  velocity,  retreating  side  of  the  rotor  disk.  The  low  velocities  require 
retreating  blades  to  be  placed  at  high  angles)  of  attack  well  into  the  stall  regime. 
Due  to  the  dynamic  nature  of  the  stall  penetration,  static  airfoil  characteristics 
have  been  found  inadequate  in  predicting  force  characteristics  which  are  experi¬ 
enced  by  stalled  rotor  blades.  To  predict  the  rotor  blade  forces  requires  utiliza¬ 
tion  of  an  unsteady  aerodynamic  theory  such  as  developed  in  this  report. 

The  following  model  and  full-scale  test  data  serve  to  demonstrate  the  degree  to 
which  the  unsteady  aerodynamics  and  yawed  flow  refinements  improved  overall 
rotor  prediction  capability  of  the  Government's  Rotorcraft  Flight  Simulation 
Program,  C-81. 

ISOLATED  ROTOR  CORRELATION 

Figures  26  through  37  present  co-relation  of  the  G-81  program  with  available  test 

25 

data  for  an  8-foot-diameter  model  wind  tunner  rotor  .  The  model  rotor  repre¬ 
sented  a  1/7.  5  scale  model  of  a  CH-47C  helicopter  rotor.  Data  are  shown  for 
an  advance  ratio  of  0. 35  and  an  advancing  tip  Mach  number  of  0. 6.  The  data  were 
obtained  during  a  rotor  shaft  angle  variation  at  constant  collective  pitch. 

Figure  26  presents  thrust  correlation  for  a  shaft  angle  sweep  at  constant  collective 
pitch.  The  theoretical  estimates  were  made  at  a  collective  pitch  angle  which 
matched  the  test  thrust  levels  of  Cij./o  =  .  074  at  the  control  axis  angle  of  -19. 8°. 
At  low  lift  levels ,  this  collective  setting  gives  thrust  predictions  which  are  below 
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ROTOR  THRU  T  COEFFICIENT  (CT/a) 
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CONTROL  AXIS  ANGLE  (dCA)  -  DEGREES 

Figure  26.  Model  Rotor  Correlation  With  Elastic  Blade  Rotor 
Theory:  Rotor  Thrust  vs.  Control  Axis  Angle. 
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Figure  28.  Model  Rotor  Correlation  With  Elastic  Blade  Rotor 
Theory:  Rotor  Thrust  vs.  Rotor  Effective  Drag. 
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Figure  30.  Model  Rotor  Correlation  With  Rigid  Blade  Rotor 

Theory:  Rotot*  Thrust  vs.  Rotor  Propulsive  Force 
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ROTOR  THRUST  COEFFICIENT  (CT/a) 
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Figure  33, 


Model  Rotor  Correlation  With  Elastic  Blade  Rotor 
Theory:  Rotor  Thrust  vs.  Alternating  Root  Torsion 
Load 
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C-81  WITH  RADIAL  FLOW  AND 
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BLADE  AZIMUTH  POSITION  -  DEGREES 

Figure  34.  Correlation  of  Theory  and  Model  Rotor  Blade  Root 
Torsion  Loads:  Blade  Root  Torsion  Load  vs.  Blade 
Azimuth  Position. 
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Figure  35.  Correlation  of  Theory  and  Model  Rotor  Blade  Root 
Torsion  Loads:  Blade  Root  Torsion  Load  vs.  Blade 
Azimuth  Position, 
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test  levels.  At  high  lift  levels,  the  opposite  condition  exists:  the  theory  overpre¬ 
dicts  thrust  levels,  This  condition  basically  results  in  a  discrepancy  in  the  ability 
of  the  C-81  program  to  predict  accurate  lift  slope  for  the  shaft  sweep.  This  type 
of  discrepancy  between  theory  and  test  has  been  found  by  other  investigators1  for 
other  rotor  analyses.  (See  for  example  correlations  presented  in  Reference  3.) 

The  only  significant  result  of  this  circumstance  is  that  the  analysis  cannot  predict 
accurate  collective  estimates  at  extremely  low  thrust  levels.  However,  later 
figures  will  illustrate  that  for  low  thrust  levels,  the  performance  estimates  are 
fairly  accurate. 

The  unsteady  aerodynamics  and  radial  Clow  methods  developed  under  this  contract 
produce  no  differences  in  the  prediction  of  thrust  versus  control  axis  at  low  thrust 
levels.  However,  at  high  thrust  levels,  the  effects  of  unsteady  aero.and  radial  flow 
provide  improved  prediction  capability.  The  basic  C-81  analysis  is  unable  to  pre¬ 
dict  the  high  thrust  levels  which  are  measured  during  ro  >r  model  tests.  The  addi¬ 
tion  of  radial  flow  pro  .ides  increased  lift  prediction  at  high  thrust  levels  due  to 
extended  Ci  capability  added  by  the  radial  flow  theory.  This  theory  is  still 
deficient  at  the  extreme  high  thrust  levels  of  the  test  data.  Adding  the  unsteady 
aerodynamics  theory  allows  the  theory  to  predict  the  test  thrust  levels  and  in 
addition  provides  improved  correlation  of  the  thrust  levels  and  lift  slope  in  the 
range  of  C>p/ o  of  .  08  and  above.  Utilization  of  the  combined  unsteady  aerodynamic 
and  radial  flow  theories  in  C-81  provides  the  best  correlation  with  test  data  above 
stall. 

Rotor  performance  data  corresponding  to  the  data  of  Figure  26  are  presented  in 
Figures  27  and  28.  Again,  the  combined  unsteady  aerodynamics  and  radial  flow 
theories  provide  the  most  accurate  prediction  of  the  test  results.  At  low  lift  levels, 
the  effects  of  increased  skin  friction  drag  levels  due  to  radial  flow  adjust  the  pre- 
dic  ed.  rotor  propulsive  force  coefficient,  X/qd2o  ,  and  rotor  effective  drag  coeffi¬ 
cient,  De/qd2o,  to  levels  obtained  during  the  rotor  test.  At  lift  levels  where  stall 
becomes  significant,  the  combined  unsteady  aerodynamics  and  radial  flow  Ineories 
provide  the  most  accurate  estimate  of  the  data  trends.  Again,  the  chief  factor  is  the 
ability  of  the  combined  theory  to  predict  more  accurate  thrust  levels  compared  to 
the  test  data. 

The  correlations  presented  in  Figures  26,  27,  and  28  are  based  on  C-81  estimates 
which  were  obtained  using  the  elastic  rotor  option  of  the  program.  Figures  29,  30 


and  31  present  the  same  correlation  with  the  exception  that  the  elastic  rotor  capa¬ 
bility  was  not  utilized,  i.  e. ,  the  blades  are  represented  as  being  rigid  and  are  free 
to  flap  about  the  flapping  hinge.  The  correlation  results  are  similar  to  those  ob¬ 
tained  in  the  previous  correlations  with  the  elastic  rotor  representation.  The  only 
significant  difference  occurs  at  the  highest  thrust  conditions  where  the  correlation 
indicates  that  the  rigid  blade  analysis  is  unable  to  predict  the  extreme  high  thrust 
levels  of  the  test  rotor.  This  indicates  that  elastic  rotor  representation  is  required 
to  obtain  accurate  prediction  of  rotor  thrust  in  deep  stall. 

Elastic  blade  torsion  response  of  the  8-foot  model  rotor  is  shown  in  Figures  32  and 
33.  Steady  components  of  the  torsional  loads  are  predicted  well  by  the  revised 
C-81  program.  The  improved  correlation  is  attributed  both  to  improving  the 
dynamic  representation  of  blade  element  pitching  moment,  and  the  ability  of  the 
combined  unsteady  aero  and  radial  flow  theories  to  predict  high  rotor  thrust 
levels.  Similar  correlation  results  are  obtained  for  the  alternating  torsion  loads 
which  are  presented  in  Figure  33.  Typical  torsion  waveforms  are  given  in  Figures 
34  and  35  for  a  case  below  rotor  stall  and  a  case  above  rotor  stall.  The  improved 
C-81  program  provides  better  correlation  in  both  examples. 

Figures  36  and  37  compare  measured  and  theoretical  Cj  distributions  for  a  rotor 
blade  element  at  .75  radius  as  a  function  of  azimuth  position.  The  8-foot  diam¬ 
eter  model  rotor25  was  equipped  with  an  instrumented  blade,  capable  of  measuring 
chordwise  normal  pressure  distributions.  The  pressure  distributions  were  then 
integrated  over  the  chord  to  obtain  section  Cn  and  Cc*  Section  Cj  was  obtained  by 
combining  Cm  and  Cc  in  the  following  equation* 

Cj  =  Cn  cos  a  ••  Cc  sin  a  (82) 

For  an  unstalled  thrust  coefficient  of  .074,  Figure  36  illustrates  little  influence 
due  to  radial  flow  and  unsteady  aerodynamics.  For  a  thrust  coefficient  of  .  114, 
Figure  37,  where  the  rotor  is  operating  in  a  stalled  condition,  the  correlation  is 
considerably  improved  with  the  use  of  the  radial  flow  and  unsteady  aerodynamics 
option  of  the  computer  program.  The  primary  influence  is  obtains  between  azi  ¬ 
muth  positions  of  200°  and  260°  where  the  improved  analysis  more  closely  predicts 
the  high  Ci  levels  achieved  during  the  rotor  test  The  excessive  Cj  variation  be¬ 
tween  azimuth  positions  of  260°  and  330°  is  due  to  excessive  torsional  response  of 
the  predicted  blade  mr  ‘on  at  high  thrust  levels,  which  was  previously  shown  in 
Figure  33. 
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CORRELATION  OF  THEORY  AND  HELICOPTER  TEST  DATA 


UH-1H  Single-Rotor  Helicopter 

Figure  38  compares  the  C-81  analysis  with  the  UH-1H  test  data.  This  correlation 
was  performed  using  fully  coupled  elastic  mode  shapes  for  the  UII-1H  rotor  repre¬ 
sentation.  As  the  figure  indicates,  the  correlation  is  excellent. for  both  the  basic 
C-81  and  the  revised  program.  Even  though  a  9500-lb  gross  weigut  is  near  the 
ma’^'num  weight  of  the  UH-1H,  this  case  has  a  thrust  coefficient  ,(Cx/°  )  of  only 
0. 08.  For  the  maximum  test  advance  ratio  of  .25,  this  thrust  level  does  not  rep¬ 
resent  a  case  in  which  rotor  stall  is  significant.  Therefore,  the  use  of  unsteady 
aerodynamics  and  radial  flow  methods  at  low  lift  levels  does  not  have  a  significant 
effect  on  the  rotor  performance  of  a  UH-1H  configuration. 

CH-47C  Tandem  Helicopter 

The  CH-47C  helicopter  is  capable  of  attaining  advance  ratios  and  gross  weights 
where  the  effects  of  rotor  stall  become  significant.  Figure  39  illustrates  typical 
correlation  at  gross  weights  of  30,000  lb  (Cy/cr  =  .  056)  and  50,000  lb  (C-p/o  = 

.  094).  These  correlations  were  performed  without  using  the  elastic  rotor  option 
of  the  C-81  program.  This  was  done  since  the  basic  Government-supplied  C-81 
analysis  was  found  to  be  numerically  unstable  when  attempting  to  analyze  a  fully 
articulated  rotor  configuration  with  elastic  mode  shapes.  The  instability  precluded 
performing  the  correlation  with  the  use  of  elastic  rotor  .characteristics.  For  the 
30,000-lb  gross  weight  case,  the  use  of  the  combined  unsteady  aerodynamics  and 
radial  flow  theories  gives  good  correlation  with  the  test  data.  At  50, 000  lb  the  re¬ 
vised  theory  provides  improved  correlation;  however,  at  speeds  above  120  knots 
the  test  values  are  still  higher  than  the  estimated  results.  For  this  case  (C-p/o  of 
.094),  the  rotor  would  experience  significant  stall  effects  and  the  elastic  response 
would  cause  additional  losses  at  high  speeds.  It  is  anticipated  that  including  an 
elastic  rotor  representation  would  provide  a  better  rotor  estimate  for  the  50,000-lb 
condition  at  high  speeds. 

Data  shown  in  Figure  39  are  for  the  forward  rotor  performance  only.  The  correla¬ 
tions  shown  here  are  indicative  of  single-rotor  configuration  or  of  a  forward-rotor 
on  a  tandem-rotor  configuration,  since  it  is  essentially  an  isolated  rotor.  Figure 
40  illustrates  that  the  C-81  analysis  does  not  consider  the  rotor-rotor  interference 
effects  which  are  encountered  by  the  aft  rotor  of  a  tandem  configuration. 
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Figure  39.  Correlation  of  Rigid  Blade  Rotor  Theory  With 
CH-47C  Force  and  Rotor  Test  Data. 


Figure  40  includes  both  the  forward  and  aft  rotor  correlations  for  the  CH-47C 
helicopter  at  a  gross  weight  of  30, 000  lb.  Since  the  C-81  analysis  does  not -include 
an  aft  rotor  interference  penalty,  a  first  approximation  correction  has  been  applied 
in  order  to  illustrate  the  resulting  effect  on  performance  level.  The  forward  rotor 
of  a  tandem  configuration  produces  an  interference  velocity  which  acts  on  the  aft 
rotor.  Due  to  the  interference  velocity,  the  aft  rotor  essentially  operates  in  a 
climb  condition.  The  approximation  used  here  combines  the  interference  velocity 
(VINt)  wit^  r°tor  thrust  to  obtain  an  interference  pow'er  correction.  The 

interference  power  correction  is  computed  by 


HP 


INT 


VINT  (TAFT  ROTOR* 
550 


(83) 


The  detailed  effects  on  the  aft  rotor  trim  are  outside  the  scope  of  this  report,  but 
should  be  considered  when  analyzing  aft  rotor  flight  conditions.  Derivation  of  the 
interference  velocity  (VjNT)  is  based  on  considerable  theoretical  analyses  and 
model  tests  conducted  by  The  Boeing  Company.  With  the  correction  applied,  the 
C-81  analysis  is  shown  to  give  improved  correlation  when  the  unsteady  aerodynam¬ 
ics  and  radial  flow  options  are  included. 

A  CH-47C  collective  pull-up  maneuver  is  illustrated  in  Figure  40.  The  maneuver 

28 

was  performed  as  part  of  the  CH-47C  structural  demonstration  program  .  From 
an  initial  level  flight  trim  condition,  collective  pitch  was  decreased  and  then 
rapidly  increased  to  100%  collective  in  order  to  demonstrate  the  "g"  capability  of 
the  aircraft.  Figure  41  compares  the  load  factor  (n)  time  history  for  the  lest 
maneuver  with  theoretical  estimates  of  C-81.  Both  versions  of  the  program  give 
adequate  correlation  with  the  test  data;  however,  the  combined  unsteady  aerody¬ 
namics  and  radial  flow  theory  is  able  to  predict  a  slightly  higher  peak  "g"  level 
than  the  basic  C-81  program.  Again,  the  elastic  rotor  representation  could  not  be 
utilized  for  the  CH-47C  configuration  due  to  numerical  instability  of  the  basic  C-81 
program.  Based  on  the  results  of  Figures  26  and  31,  which  show  that  a  substantial 
increase  in  maximum  thrust  is  predicted  when  using  elastic  rotor  characteristics, 
it  is  anticipated  that  the  revised  C-81  program  would  predict  further  increases  in 
"g"  level  if  elastic  rotor  characteristics  could  be  used.  The  oscillation  of  the 
estimated  values  corresponds  to  a  3/rev  variation  of  rotor  thrust.  This  vibratory 
frequency  of  a  three-bladed  rotor  does  not  appear  in  the  test  data  because  it  was 
electronically  filtered  out  by  the  test  instrumentation. 
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APPENDIX  m 

SAMPLE  INPUT  AND  OUTPUT 


Input  and  output  formats  of  the  revised  C-81  program  are  similar  to  the  C-Sl 
(AGAJ7 1) ,  Reference  24.  Changes  in  input  requirements  consist  of  six  additional 
inputs  for  the  forward  rotor  and  aft  rotor  characteristics.  The  precise  definitions 
of  these  additional  inputs  are  described  in  the  USER'S  MANUAL  section.  Figure 
42  illustrates  a  typical  input  listing.  The  new  inputs  are  located  immediately  after 
the  MAIN  ROTOR  GROUP  and  immediately  after  the  TAIL  ROTOR  GROUP. 
Headings  for  the  new  inputs  are: 

YAW  FLOW:  on/off  switch  for  radial  flow  effects. 

ALFACL=0(1)  through  ALFACL=0(4):  coefficients  for  a  third-order  curve  fit 
which  represents  section  angle  of  attack  at  zero  lift  coefficient  as  a  function  of 
Mach  number. 

T/C:  thickness  ratio  of  the  rotor  airfoil  section. 

(For  detailed  definition  of  all  other  inputs,  consult  Reference  24.) 

All  output  formats  for  rotorcraft  trim  computations  remain  unchanged  from  those 
described  in  Reference  24.  Figure  43  presents  a  typical  final  trim  output  listing. 

For  maneuver  cases,  the  output  also  remains  unchanged  with  the  exception  of  the 
page  for  Rotor  Airload  Data.  Figure  44  presents  a  typical  airloads  output.  The 
revised  output  consists  of: 

PSI:  blade  azimuth  position 

U-HUB,  V-HUB,  W-HUB:  hub  velocity  components 
GEO.  PITCH:  blade  geometric  pitch 
BETA  (HUB):  blade  flap  angle  referenced  to  hub 
STA:  blade  station 

UT,  UP,  UR:  blade  element  local  velocity  components 

MACH:  blade  element  local  Mach  number 

ALPHA:  blade  element  angle  of  attack 

ALPHAD:  time  derivative  of  blade  element  angle  of  attack 


125 


THETAD:  time  derivative  of  blade  pitch  angle 

THETADD:  second  time  derivative  of  blade  pitch  angle 

CL,  CD,  CM:  blade  element  lift,  drag  and  moment  coefficients  (including 
unsteady  aero  and  radial  flow  effects ,  if  used) 

CDSF:  skin  friction  drag  coefficient  used  in  radial  flow  computations 
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